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Q . Abstract 

' We investigate the singular sets of solutions of a natural family of 

conformally covariant pscudodiffcrcntial elliptic operators of fractional 
order, with the goal of developing generalizations of some well-known 
properties of solutions of the singular Yamabe problem. 

<N ' 

^ ! 1 Introduction 

m . 

. Let (M,g) be a compact ra-dimensional Riemannian manifold, n > 3. If 

A C M is any closed set, then the 'standard' singular Yamabe problem 
concerns the existence and geometric properties of complete metrics of the 

4 

form g = u n ~ 2 g with constant scalar curvature. This corresponds to solving 
. the partial differential equation 

><: ^ + i^ ) R * u = w 1 T ) R ' um ' "><>• (L1) 

where R 9 is constant and with a 'boundary condition' that u — > oo suffi- 
ciently quickly at A so that g is complete. (Note that in our convention, A^- 
is an operator with nonnegative spectrum.) It is known that solutions with 
R 9 < exist quite generally if A is large in a capacitary sense |17| . whereas 
for R 9 > existence is only known when A is a smooth submanifold (possibly 
with boundary) of dimension k < (n — 2)/2, see [19], [8]. 

There are both analytic and geometric motivations for studying this 
problem. For example, in the positive case (R 9 > 0), solutions to this 
problem are actually weak solutions across the singular set, so these results 
fit into the broader investigation of possible singular sets of weak solutions of 
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semilinear elliptic equations. On the geometric side is a well-known theorem 
by Schoen and Yau [25] stating that if (M, h) is a compact manifold with 
a locally conformally flat metric h of positive scalar curvature, then the 
developing map D from the universal cover M to S n , which by definition 
is conformal, is injective, and moreover, A := S n \ D(M) has Hausdorff 
dimension less than or equal to (n— 2)/2. Regarding the lifted metric h on M 
as a metric on f2, this provides an interesting class of solutions of the singular 
Yamabe problem which are periodic with respect to a Kleinian group, and 
for which the singular set A is typically nonrectifiable. More generally, that 
paper also shows that if g is the standard round metric on the sphere and if 

4 _ 

g = u n ~ 2 g is a complete metric with positive scalar curvature and bounded 
Ricci curvature on a domain Q = S n \ A, then dim A < (n — 2)/2. 

In the past two decades it has been realized that the conformal Lapla- 
cian, which is the operator appearing as the linear part of (jl.ip . fits into a 
holomorphic family of conformally covariant elliptic pseudodifferential op- 
erators. The operators in this family of positive even integer order are the 
GJMS operators, and these have a central role in conformal geometry. Just 
as the Yamabe problem is naturally associated to the conformal Laplacian, 
so too are there higher order Yamabe-type problems associated to the other 
GJMS operators, or more generally, also to these other conformally covari- 
ant operators with noninteger order. The higher (integer) order Yamabe 
problems have proved to be analytically challenging and provide insight 
into the GJMS operators themselves. Hence it is reasonable to hope that 
these fractional order (singular) Yamabe problems will have a similarly rich 
development and will bring out interesting features of these conformally co- 
variant pseudodifferential operators. Prom a purely analytic perspective, 
little is known about regularity of solutions of semilinear pseudodifferential 
equations like these, and this family of geometric problems is a natural place 
to start. In fact, as we explain below, fractional powers of the Laplacian 
have also appeared recently in the work of Caffarelli and his collaborators as 
generalized Dirichlet to Neumann operators for certain singular divergence 
form elliptic equations, which further indicates the worth of studying such 
operators. 

The present paper begins an investigation into these questions. Our 
goals here are limited: beyond presenting this set of problems as an in- 
teresting area of investigation, we prove a few results which indicate how 
certain properties of the fractional singular Yamabe problem extend some 
well-known results for the standard Yamabe equation. 

To describe this more carefully, we first define the family of fractional 
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conformal powers of the Laplacian. As we have already indicated, the linear 
operator which appears as the first two terms on the left in (jl.ip is known 
as the conformal Laplacian associated to the metric ~g, and denoted Pf . It 
is conformally covariant in the sense that if / is any (smooth) function and 

4 _ 

g = u n ~ 2 g for some u > 0, then 

P?(uf) = u^Pf(f). (1.2) 

Setting / = 1 in (jl.2p yields the familiar relationship (jl.ip between the 
scalar curvatures R 9 and RP . P\ is the first in a sequence of conformally 
covariant elliptic operators, Pk, which exist for all k £ N if n is odd, but only 
for k G {1, . . . ,n/2} if k is even. The first construction of these operators, 
by Graham-Jenne-Mason-Sparling [13] (for which reason they are known as 
the GJMS operators), proceeded by trying to find lower order geometric 
correction terms to A fc in order to obtain nice transformation properties 
under conformal changes of metric. Beyond the case k = 1 which we have 
already discussed, the operator 

P 2 = A 2 + 5 (a n Rg + b n Ric) d + ^Q 2 , 

called the Paneitz operator (here Q2 is the standard Q-curvature), had also 
been discovered much earlier than the operators P^ with k > 2. 

This leads naturally to the question whether there exist any conformally 
covariant pseudo differential operators of noninteger order. A partial result 
in this direction was given by Peterson [20], who showed that for any 7, the 
conformal covariance condition determines the full Riemannian symbol of a 
pseudodifferential operator with principal symbol |£| 27 . Hence P 7 is deter- 
mined modulo smoothing operators, but it is by no means clear that one can 
choose smoothing operators to make the conformal covariance relationships 
hold exactly. The breakthrough result, by Graham and Zworski [T3], was 
that if (M, [g]) is a smooth compact manifold endowed with a conformal 
structure, then the operators Pk can be realized as residues at the values 
7 = k of the meromorphic family S(n/2 + 7) of scattering operators associ- 
ated to the Laplacian on any Poincare-Einstein manifold (X, G) for which 
(M, [g]) is the conformal infinity. These are the 'trivial' poles of the scatter- 
ing operator, so-called because their location is independent of the interior 
geometry; S(s) typically has infinitely many other poles, which are called 
resonances, the location and asymptotic distribution of which is a matter 
of considerable interest and ongoing study. Multiplying this scattering fam- 
ily by some T factors to regularize these poles, one obtains a holomorphic 
family of elliptic pseudodifferential operators P^ (which patently depends 
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on the filling (X, G)). An alternate construction of these operators has been 
obtained by Juhl, and his monograph [15] describes an intriguing general 
framework for studying conformally covariant operators, see also [16] . 

This realization of the GJMS operators has led to important new under- 
standing of them, including for example the basic fact that P® is symmetric 
with respect to dVj, (something not obvious from the previous fundamen- 
tally algebraic construction). Hence even though the family P 7 is not en- 
tirely canonically associated to (M, [g]) (as we explain in some detail below), 
its study can still illuminate the truly canonical operators which occur as 
special values at positive integers, i.e. the GJMS operators. 

For various technical reasons, we focus here only on the operators P 7 
when 7 6 R, I7I < n/2. These have the following properties: first, Pq = Id, 
and more generally, Pk is the /c th GJMS operator, k = 1, . . . , n/2; next, P 7 
is a classical elliptic pseudodifferential operator of order 2-y with principal 
symbol G2 1 {P^) = ', hence (since M is compact), P 7 is Fredholm on L 2 
when 7 > 0; if P 7 is invertible, then P_ 7 = PZ~ ; finally, 

4 _ n + 2~i 

if g = u —g, then P°(uf) = u~P°(f) (1.3) 

for any smooth function /. Generalizing the formulas for scalar curvature 
(7 = 1) and the Paneitz-Branson Q-curvature (7 = 2), we make the defini- 
tion that for any < 7 < n/2, Q 7 , the Q-curvature of order 7 associated to 
a metric g, is given by 

Q? = -P?(l)- (1-4) 
Let us comment further on the choices involved in these definitions. 
First, Poincare-Einstein fillings (X,G) of (M, [g]) (which are defined at the 
beginning of §2), may not always exist, and when they exist, they may not 
be unique. The existence issue is not serious: the construction of [H] only 
uses that the metric G satisfy the Einstein equation to sufficiently high or- 
der, and one can even take X = M x [0, 1] with the conformal structure \g\ 
at M x {0} and with the other boundary M x {1} a regular (incomplete) 
boundary for G. However, these comments indicate that the issue of lack of 
uniqueness is far worse, since there are always infinite dimensional families 
of asymptotically Poincare-Einstein fillings. Any choice of one of these fixes 
a family of operators P 7 , and for each such choice P 7 satisfies all the proper- 
ties listed above. As already noted, the complete Riemannian symbol of P 7 
is determined by the metric ~g and the conformal covariance; the choice of 
filling provides a consistent selection of smoothing terms in these pseudod- 
ifferential operators for which the same covariance properties hold. Hence 
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the Q-curvatures Q 7 for noninteger values of 7 are similarly ill-defined. In 
particular, except in certain special cases where there are canonical choices 
of fillings (e.g. the sphere), it is not clear that the existence of a metric ~g in a 
conformal class such that Q 7 > depends only on that conformal class. We 
leave open these significant problems, and in what follows, always make the 
tacit assumption that for any given (M, [<?]), we have fixed an approximately 
Poincare-Einstein filling (X,G) and used this to define the family P 7 . In 
other words, it is perhaps more sensible to think of P 7 and Q 7 as quantities 
determined by the pair ((M, [g]), (X, G)). 

In any case, generalizing (jl.ip . consider the "fractional Yamabe prob- 
lem": given a metric ~g on a compact manifold M, find u > so that if 
g = ti 4// ( n_27 ^, then Q 7 is constant. This amounts to solving 

_ 11+27 
FS,u = Q»u»-*r , u>0, (1.5) 

for Q 7 = const. More generally, we can simply seek metrics g which are 
conformally related to g and such that Q 7 > or Q 7 < everywhere. 

This fractional Yamabe problem has now been solved in many cases 
where the positive mass theorem is not needed [12] . and further work on 
this is in progress. 

As described earlier, it is is also interesting to construct complete metrics 
of constant (positive) Q 7 curvature on open subdomains Q = M \ A, or in 
other words, to find metrics g = u^^^^g which are complete on and 
such that u satisfies (|1.5p with Q 7 a constant. This is the fractional singular 
Yamabe problem. In the first few integer cases it is known that the positivity 
of the curvature places restrictions on dim A (see (25], [19] for the case 7 = 1, 
[3] for 7 = 2, and [11] for the analogous problem for the closely related 
curvature) . 

Although it is not at all clear how to define P 7 and Q 7 on a general com- 
plete open manifold, we can give a reasonable definition when f2 is an open 
dense set in a compact manifold M and the metric g is conformally related 
to a smooth metric g on M. Namely, we can define them by demanding that 
the relationship (ll.3p holds. Note, however, that this too is not as simple 
as it first appears since, because of the nonlocal character of P 7 , we must 
extend u as a distribution on all of M. We discuss this further below. 

The purpose of this note is to clarify some basic features of this fractional 
singular Yamabe problem and to establish a few preliminary results about 
it. Our first result generalizes the Schoen-Yau theorem. 

_ 4 
Theorem 1.1. Suppose that (M n ,g) is compact and g = u n - 2 ig is a com- 
plete metric on £1 = M\A, where A is a smooth k-dimensional submanifold. 
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Assume furthermore that u is polyhomogeneous along A with leading expo- 
nent —n/2 + 7. If < 7 < 7^, and if > everywhere for any choice 
of asymptotically Poincare- Einstein extension (X, G) which defines P% and 
hence Q^, then n, k and 7 are restricted by the inequality 

r (i-M)A(i-H)>°. ™ 

where T is the ordinary Gamma function. This inequality holds in particular 
when k < (n — 2^)/2, and in this case then there is a unique distributional 
extension of u on all of M which is still a solution of (|1.5p . 

Remark 1. Recall that u is said to be polyhomogeneous along A if in terms 
of any cylindrical coordinate system (r,9,y) in a tubular neighborhood of A, 
where r and 6 are polar coordinates in disks in the normal bundle and y is 
a local coordinate along A, u admits an asymptotic expansion 

where fij is a sequence of complex numbers with real part tending to infinity, 
for each j, is nonzero for only finitely many nonnegative integers k, and 
such that every coefficient ajk € C°° . The number fiQ is called the leading 
exponent if > K(/Uo) for all j ^ 0. We refer to fW/ for a more 

thorough account of poly homogeneity. 

Remark 2. As we have noted, inequality (jl.6j) is satisfied whenever k < 
(n — 2j)/2, and in fact is equivalent to this simpler inequality when 7 = 1. 
When 7 = 2, i.e. for the standard Q— curvature, this result is already known: 
it is shown in ^ that complete metrics with Q2 > and positive scalar 
curvature must have singular set with dimension less than (n — 4)/2, which 
again agrees with (jl.6p , 

We also present a few special existence results. First, the following re- 
mark exhibits solutions coming from Kleinian group theory where A is non- 
rectifiable. 

Remark 3. Suppose that 7 G [1, n/2). Let T be a convex cocompact subgroup 
of SO(n + 1, 1) with Poincare exponent 5(F) G [1, (n - 27)/2). Let A C S n 
be the limit set ofT. Then 17 = S n \A admits a complete metric g conformal 
to the round metric and with > 0. 

As we explain below, this follows directly from the work of Qing and 
Raske [22]. 

Finally, one can also obtain existence of solutions when 7 is sufficiently 
near 1 and A is smooth by perturbation theory. 
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Theorem 1.2. Let (M n , [g]) be compact with nonnegative Yamabe constant 
and A a k- dimensional submanifold with k < ^(n — 2). Then there exists an 
e > such that 1/7 G (1 — e, 1 + e), there exists a solution to the fractional 
singular Yamabe problem (11.51) with > which is complete on M \ A. 

Our final result is a growth estimate for weak solutions that are singular 
on S n \£l. Our result is not very strong in the sense that we do need to 
require that u is a weak solution in the whole S n . However, it provides the 
first insight into a general theory of weak solutions on subdomains of S n . 

Proposition 1.3. Let g c be the standard round metric on S n , and (B n+1 , G) 

the Poincare ball model of hyperbolic space, which has (S n , [g c ]) as its con- 

4 

formal infinity. Let g = u n ~ 2 i g c be a complete metric on a dense subdomain 
of the sphere, SI = S n \ A, with Q 7 equal to a positive constant, and such 
that u is a distributional solution to 

n + 2~f 

P9c u = U"~ 2 ~' (1.7) 

on S n (with u finite only ontl). Then, for all z £ Q, 

C 



u(z) < 



, / an "~ 2 7 

dg c (z,A) 2 



where C depends only on n and 7. 

There are many interesting questions not addressed here. For example, 
we point out again that there is not yet a good definition of the family P® 
on an arbitrary complete manifold (fi, g). Provided one is able to make this 
definition, it would then be useful to compute the L 2 -spectrum of P%, even 
for some specific examples such as HP or W k+l x S n ~ k ~ 1 . Finally, it would 
also be important to obtain the correct generalization of the Schoen-Yau 
theorem for the operators P 7 . We hope to address these and other problems 
elsewhere. 



2 Fractional conformal Laplacians 

We now provide a more careful description of the construction of the family 
of conformally covariant operators P 7 , and also give two alternate definitions 
of these operators in the flat case to provide some perspective. 

As we have described in the introduction, Graham and Zworski p3] dis- 
covered a beautiful connection between the scattering theory of the Lapla- 
cian on an asymptotically hyperbolic Einstein manifold and the GJMS op- 
erators on its conformal infinity. Let (M,g) be a compact n-dimensional 
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Riemannian manifold. Suppose that X is a smooth compact manifold with 
boundary, with dX = M, and denote by re a defining function for the bound- 
ary, i.e. x > on X, x = precisely on dX and dx 7^ there. A metric G 
on the interior of X is called conformally compact if x 2 G = G extends as a 
smooth nondegenerate metric on the closed manifold with boundary. It is 
not hard to check that G is complete and, provided that = 1 at dX, 

the sectional curvatures of G all tend to —1 at 'infinity'. The metric G is 
called Poincare-Einstein if it is conformally compact and also satisfies the 
Einstein equation Ric G = —nG. As we have explained, it is only necessary 
to consider asymptotically Poincare-Einstein metrics; by definition, these 
are conformally compact metrics which satisfy Ric G = —nG + G(x N ) for 
some suitably large N (typically, N > n is sufficient). 

The conformal infinity of G is the conformal class of G\ TQX ; only the 
conformal class is well defined since the defining function x is defined up to 
a positive smooth multiple. If g is any representative of this conformal class, 
then there is a unique defining function x for M such that G = x~ 2 (dx 2 + 
g(x)) where g(x) is a family of metrics on M (or rather, the level sets of x), 
with g(0) the given initial metric. 

We now define the scattering operator S(s) for (X,G). Fix any /o € 
C°°(M); then for all but a discrete set of values s£C, there exists a unique 
generalized eigenfunction u of the Laplace operator on X with eigenvalue 
s(n — s). In other words, u satisfies 

( (A G - s(n - s))u = 

\ u = f x ns + /V, for some L f e C oo (X) with / | a=o = /q 

By definition, S(s)fo = f\ x =o- This is an elliptic pseudodifferential operator 
of order 2s — n which depends meromorphically on s; it is known to always 
have simple poles at the values s = n/2, n/2 + 1, n/2 + 2, . . .. These locations 
are independent of (X, G), hence are called the trivial poles of the scattering 
operator. S(s) has infinitely many other poles which are of great interest in 
other investigations, but do not concern us here. Letting s = n/2 + 7, we 
now define 

J »_J»rm s (» +7 ) ! (2.9) 
because of these prefactors, one has that the principal symbol is 

^2 7 (P|) = \v\p- (2-10) 

The scattering operator satisfies a functional equation, S(s)S(n — s) = Id, 
which implies that 

P 7 oP„ 7 = Id. (2.11) 
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Finally, it is proved in |14j that the operators P 7 satisfy the conformal 
covariance equation ([1.3D , 

This definition of the operators P 7 depends crucially on the choice of the 
Poincare-Einstein filling (X,G). Graham and Zworski point out that it is 
only necessary that the metric G satisfy the Einstein equation to sufficiently 
high order as x — >• in order that the properties of the P 7 listed above be 
true (for 7 in a finite range which depends on the order to which G satisfies 
the Einstein equation). As we have discussed in the introduction, it is always 
possible to find such metrics, and we suppose that one has been fixed. 

Let us now address the issue of how to define P 7 and Q 7 when £1 is a 
dense open set in a compact manifold M and g is complete and conformal 
to a metric g which extends to all of M. (As usual, we assume that (M, ~g) 
has an asymptotically Poincare-Einstein filling). There is no difficulty in 
using the relationship (jl.3p to define Pljf when / € C^°(0) . From here one 
can use an abstract functional analytic argument to extend P 7 to act on 
any / € L 2 {£l,dV g ). Indeed, it is straightforward to check that the operator 
P 7 defined in this way is essentially self-adjoint on L 2 (£l,dV g ) when 7 is 
real. However, observe that P 7 = A^ + K, where K is a pseudo-differential 
operator of order 27 — 1. Furthermore, A 7 is self-adjoint by the functional 
calculus, so we can appeal to a classical theorem, see [23], which states that 
a lower order symmetric perturbation of a self-adjoint operator is essentially 
self-adjoint. 

A separate, but also very interesting issue, is whether Q 7 is a positive 
constant implies that the conformal factor u is a weak solution of (|1.5|) on 
all of M. This is true (with some additional hypotheses) when 7 = 1, cf. 
[25]. 

We conclude this section with two alternate definitions of the operators 
P 7 in the special case where (M, [g]) = 1" with its standard flat conformal 
class. 

The canonical Poincare-Einstein filling in this case is the hyperbolic space 
X = = R+ x M™ with metric G = x~ 2 (dx 2 + \dy\ 2 ). 

Since g is flat, we have P 7 = Aj, and this can be written in either of the 
two equivalent forms: 

Ay(y) = (27T)-™ / e^\ V \ 2 y( V )d V or 

f(y) - f(y) 
\y-y\ 



Both formula? can be regularized so as to hold for any given 7. 
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One other way that A 7 arises is as a generalized Dirichlet to Neumann 
map; this definition is essentially the same as the one above involving the 
scattering operator (indeed, the point of view in geometric scattering theory 
is that the scattering operator is simply the Dirichlet to Neumann operator 
at infinity), but as recently rediscovered by Chang-Gonzalez [3] in relation to 
the work on (Euclidean) fractional Laplacians by Caffarelli and Silvestre [2] , 
it is sometimes helpful to consider the equation in a slightly different form. In 
the following result, let (X, G) be an asymptotically Poincare-Einstein filling 
of the compact manifold (M n , [g]). Fix a representative g of the conformal 
class on the boundary and let x be the boundary defining function on X 
such that g = x~ 2 (dx 2 + g x ) with go = g. Also, write G = x 2 G; this is an 
incomplete metric on X which is smooth (or at least polyhomogeneous) up 
to the boundary. 

Proposition 2.1. ('Sj) Let U = X2~ y u and 



Then, for any /q £ C°°(M), the eigenvalue problem (|2.8p is equivalent to 



where the divergence and gradient are taken with respect to G. Moreover, 



for some nonzero constant c? 7 depending only on 7 and n. 

The Euclidean version of this result (where (X, G) is the hyperbolic up- 
per half-space) was the one studied by Caffarelli and Silvestre. The main 
advantage in this reformulation is that certain estimates are more transpar- 
ent from this point of view. 



We now turn to the proofs of Theorems 11.11 and 11.21 and Remark [3J 





on {X, G) 
fo on M, 



(2.12) 



P|(/ ) = d 1 lim lx^W/ 



3 Proofs 
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3.1 Dimensional restrictions on singular sets 

The idea for the proof of Theorem II. II is straightforward: let u be a polyho- 
mogeneous distribution on M with singular set along the smooth submani- 
fold A. Suppose that the leading term in the expansion of u is a(y)r~ n / 2+ ~ 1 . 
Then by a standard result in microlocal analysis [6], the function P^u is 
again poly homogeneous and has leading term b(y)r~( 2 " /+n ^ 2 , where b(y) = 
Xa(y) for some constant A. Now, if u is a conformal factor for which 
g = u 4 /( n - 2 ~t)g has > 0, then Pfu > 0, which implies that A > 0. 
So we must compute A to obtain (jl.6p . 

This microlocal argument states that if u is poly homogeneous, then the 
leading term of P%u can be computed using the symbol calculus (for pseu- 
dodifferential operators and for poly homogeneous distributions), and more 
specifically, that the principal symbol of P®u is equal to the product of the 
principal symbols of P® and that of u. (Note that the principal symbol of 
a distribution conormal to a submanifold A is computed in terms of the 
Fourier transform in the fibres of NA.) In the present setting, this implies 
that the constant A is the same as for the model case when M = S n and A 
is an equatorial S k , so we now focus on this special case. 

Transform S n to R n by stereographic projection, so that A is mapped to 
a linear subspace ]R fe and g is the flat Euclidean metric (which we henceforth 
omit from the notation). Write R n 3 y = (y',y") € R k x R n ~ k , so that (in 
this model case) u(y) = \y"\~ n / 2+ ~t for the singular metric u n ~ 2 ~< g; then 

P lU {y) = A~<u(y) = (2vr)^ n f e^-^^ly"^ 2 ^ dydr] 



= (2ir) k - n / e i( - y "-v"> r >"\y"\- n / 2 +7dy"d< n ". 
JM. n - k xR n - k 

Now recall a well-known formula for the Fourier transform of homoge- 
neous distributions in R^: 



e 



where 

c(N a) = vr""^ 2 F(a/2) 
{ ' J r((JV-a)/2) 

Applying this formula with N = n — k (and replacing y" by y and rj" by 77, 
for simplicity) yields first that 

e -^| y |-"/2+7 dy = c ( n _ ki !l_ k + 1 \ |„|"f +*-T 
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then, multiplying by \n\ 2 ^ and taking inverse Fourier transform we obtain 

1 / n \ ( n \ , ,_n 

- — r — r c[n — k, k + 7 c ra — k, — h 7 u 2 ~. 

(27r) n - fc V '2 J V ' 2 / 11 

Altogether then, the multiplicative factor A is equal to 

r(|(f-7)) r(£(f-fc- 7 ))- 

Discarding the factors which are always positive (which includes T(n/A— 7/2) 
since 7 < n/2), we obtain (|1.6p . 

It is unfortunately slightly messy to write down the entire set of values 
of k and 7 for which <\1.6\i holds. However, if k < (n — 2 7 )/2, then both 
§ — k ± 57 > 0. Furthermore, if 7 = 1 and A := \ - k < \, then the T 
function always takes on values with different signs at A + 7/2 and .A — 7/2. 
More generally, if we fix n and A; and let 7 increase from to n/2, then 
r(A + 7/2)/r(A — 7/2) = 1 when 7 = 0; T(A — 7/2) changes sign every time 
7 increases by 2, whereas T(j4 + 7/2) also changes similarly, but only for 7 
in the range (0, —2A), where A < 0. 

To prove the final statement of the theorem, note that if —7 — n/2, the 
leading exponent of PJju, is greater than k — n, the codimension of A, then 
-PyU cannot have any mass supported on A, which means that u is a weak 
solution of PJ/U = Q^( n+2 T)/( n - 2 T) on all of M. 



3.2 Kleinian groups 

We now turn to a special case where this problem has a direct relationship to 
hyperbolic geometry. Let T be a convex cocompact group of motions acting 
on H n+1 . Thus T acts discretely and properly discontinuously on hyperbolic 
space, is geometrically finite and contains no parabolic elements. Its domain 
of discontinuity is the maximal open set !1 C 5" on which the action of T 
extends to a discrete and properly discontinuous action; by definition of 
convex cocompactness, the quotient £l/T = Y is a compact manifold with a 
locally conformally flat structure. The complement S n \ = A is the limit 
set of r. Furthermore, the manifold X = EI n+1 /r with its hyperbolic metric 
is Poincare-Einstein with conformal infinity Y with the conformal structure 
induced from S n . We use these canonical fillings to define P 7 and Q 7 . 

In [22], Qing and Raske attack the problem of finding metrics of constant 
Qk curvature (with k < n/2 an integer). Their method involves finding met- 
rics of constant Q 7 curvature for all 1 < 7 < k. They rephrase the problem 
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P y u = Q 1 u { - n+2 ^/^ n - 2 ^ in the equivalent form u = P^(Q^ n+2 ^^ n - 2 ^). 
The advantage of this modification is that P_ 7 is a pseudodifferential oper- 
ator of negative order, and its Schwartz kernel can be obtained by summing 
the translates of the Schwartz kernel of P_ 7 on S n over the group T. This 
sum converges provided the Poincare exponent of T is less than n ~ 27 , and 
because this is a convergent sum one directly obtains explicit control and 
positivity of this operator. Fixing 1 < 7 < n/2 and restricting to convex 
cocompact groups T with Poincare exponent in this range, they are able to 
prove that if Y = fi/T has positive Yamabe type, then it admits a metric of 
constant positive Q 7 curvature. 

The proof of Remark [3] follows directly from this by lifting the conformal 
factor and solution metric to !1 C S". Namely, by the theorem of Schoen 
and Yau, the developing map of Y is injective from the universal cover Y 
to fi, and the solution metric g on Y lifts to a complete metric g on £1 of 
the form u^ n+2 "'^^ n ~ 2 " , ^g, where ~g is the standard round metric. Using the 
bound on the Poincare exponent and the compactness of Y, standard lattice 
point counting arguments show that u(p) < cdist^ (p, A)( 27_n )/ 2 . This shows 
that not only is u a solution of the modified integral equation, but is also 
a weak solution of (jl.5p on all of S n and that Q 7 is constant. Finally, by 
Patterson-Sullivan theory, the dimension of the limit set A is precisely the 
Poincare exponent 5(T). In other words, we have produced a solution to the 
fractional singular Yamabe problem with exponent 7 € [l,n/2) and with 
singular set of dimension less than n/2 — 7. 

The Qing-Raske theorem is not stated for the remaining cases 7 € (0, 1); 
it is plausible that their proof may be adapted to work then, and hence the 
lifted solution would also give a solution to our problem also for 7 in this 
range, but we do not claim this. Note, however, the results in §4 below 
concerning growth estimates for solutions of this equation for this range of 

7- 

3.3 Perturbation methods 

We come at last to the perturbation result. We deduce existence of solutions 
for the fractional singular Yamabe problem for values of 7 near 1 from the the 
general existence result in [19] for the singular Yamabe problem with 7 = 1. 
Let (M, ~g) and A be a submanifold of dimension k as in the statement of the 
theorem. (Slightly more generally, we could let the different components of A 
have different dimensions, but for simplicity we assume that A is connected.) 
Then there is a function wonM\A such that g = u 4// ( n_2 )g is complete and 
its scalar curvature Qf is a positive constant. Moreover, it is known that the 
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linearization of the equation (jl.ip at any one of the solutions u constructed 
in [19] is surjective on appropriate weighted Holder spaces. 

In the following we phrase this rigorously and parlay this surjectivity 
into an existence theorem for 7 near 1 using the implicit function theorem. 
Let T CT A denote the tube of radius a (with respect to g) around A; this 
is canonically diffeomorphic to the neighourhood of radius a in the normal 
bundle N A for a sufficiently small, and we use this to transfer cylindrical 
coordinates (r, y, 6) E [0, a) xl4 y x S n ~ k ~ 1 in a local trivialization of NA to 
Fermi coordinates in T a A. 

We use these coordinates to define weighted Holder spaces with a certain 
dilation covariance property. For w E C°(T CT A), let 

(r + r) a \w(z) — w(z)\ 

\\w\\e,0,a,0 = SUp \w\ + SUp ; ; ; — ; =— . 

zeT"A z,z&T"K \r - f\ a + \y - y\ a + (r + r) a \6 - 0\ a 

and denote by Ce'°(M\A) the space of all functions w E C°(T fT A) such that 
this norm is finite. The initial subscript e in the norm signifies that these are 
'edge' Holder spaces. Next, Ce' a {M\A) denotes the subspace of C k (M\A) 
on which the norm 

k 

\\w\\k,a,0 = \\w\\k,a,M*/2 + ^ ' ' W ' ' <M) ,q 

3=0 

is finite, where M a/2 = M\T^ T Finally, for i/Gl, let 

C^ a (M\A) = {u; = r u w : w E C e fc ' Q (M\A)} , 

with corresponding norm || • || e .fc iQ ,.^. 

Fixing Qty = 1, the linearization of u H Pfu — u ^ n+2 ^^ n ~ 2 ^ is the 
operator 

jj n + 27 4 -r 

v i—)- L y v := PJ*v u n ~ 2 ^v. 

1 n — 27 

Let u be one of the solutions to the singular Yamabe problem (7 = 1) on 
M \ A constructed in [TO]. It is proved there that the solution u has the 
form u = cir 1-n / 2 (l + v), where v E Ct ,a for any < v < k/2 and c\ > 
depends only on the dimensions k and n; furthermore, the mapping 

L 1 :C 2 u ' a (M\A)^C°^ 2 (M\A) 

is surjective for v in this same range. 
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We claim that for 7 sufficiently close to 1, and for v € (17, k/2 — rj), where 
rj > is some small fixed number, the mapping 

L 7 : C^(M \ A) — ► C° v % 2 ^ (M \ A) 

is also bounded and surjective. The boundedness follows by an interpolation 
argument. Indeed, the spaces Cq' q have interpolation properties which are 
identical to those for the ordinary Holder spaces since they are just the stan- 
dard Holder spaces for the complete metric g = g/r 2 ; a minor adjustment 
shows that the addition of the weight factor behaves as expected. The as- 
sertion about the boundedness of L 7 is clearly true for 7 = 0, 1, 2, and hence 
by interpolation is true for all 7 close to 1. (It is true for the full range of 
7 £ (0, 2) if one makes the standard change, replacing the Holder space by a 
Zygmund space, when a + 2(1 — 7) is an integer.) This also follows from |18j 
because r 27 L 7 is a pseudodifferential edge operator of order 27. Similarly, 
surjectivity follows from the construction of a parametrix for L 7 in the edge 
calculus, from [18] again. This proves that L 7 is Fredholm, and since it is 
surjective at 7 = 1, it must remain surjective for values of 7 which are close 
to 1. We write its right inverse as G 7 . 
Now consider the mapping 

(7, c, v) 1 — y JV( 7 , c, v) := G 7 (P 7 T cr 7 - n / 2 (l + v) - (cr 7 " n / 2 (l 

If U\ = C\ 

r l-n/2(l 

+vi) is the solution to the singular Yamabe problem from 
[19] , then N(l,ci,v\) = 0. Let c € (ci — s, c\ + e), and similarly, v — v\ lie 
in a ball of radius e about in Cu ■ Clearly D V N\^ L ^^ V _^ = GxLi = Id. 
The implicit function theorem now applies to show that for every (7, c) near 
to (l,ci), there exists a unique i> 7 € Ct ,a with norm less than e such that 
u 7 = cr 7 ~ n / 2 (l + w 7 ) is a solution of the fractional singular Yamabe problem 
with singular set A. 

4 Growth estimates for weak solutions on S n 

In this final section we furnish the proof of Proposition 11.31 if 7 € (0, 1) and 
H C 5" is dense, then any weak solution of the fractional singular Yamabe 
problem 

P^{u) = u"- 2 ~f in S n , u > 0, u singular along S n \fl (4.13) 

satisfies a general growth estimate. This is a direct adaptation of Schoen's 
proof (which is written out in full in [21]) for the case 7 = 1. 
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We first comment on the local regularity for solutions of (jl.7p . There are 
several ways to deduce the necessary estimates. The path we follow uses the 
equivalence, as described in l2.1| of (jl.7p with the extension problem (|2.12p : 

f -div (x^^VU) + E(x)U = 0, in(X,G), 

I -y^dJJ = Cn^U^i , on x = 0; 

here U = x^~ 7 u, G = x 2 G. 

From this point of view, we can use the linear regularity theorem [181 
Theorem 7.14] to prove that U is smooth up to x = away from A. This can 
also be deduced using standard elliptic estimates for the pseudodifferential 
operator P® c , but we refer also to more classical sources from which this can 
also be deduced, in particular the paper by [TJ; we also refer to more recent 
references [1] (where many properties of the solution are written down), and 
|12| (which holds for more general ambient metrics). In particular, from 
these last papers, one has that Schauder and local L p — > L°° estimates hold, 
and the equation also satisfies the standard maximum principles. 

Fix zq A and choose o < dist 9c (zo, A). For simplicity, write p(z) := 
distg c (z, zq). Now define 

n — 27 

f(z) := (a-p(z))— u{z); 

note that / = on dB a (xo). 

It suffices to show that f(z) < c for some c > and for all z £ B <t {zq) 
since if we choose a = dist(zo, A)/2, then f(zo) = a^^u^o), and hence 

c 

u(Zq) < 



n-27 > 

d(z ,A) 2 



which would finish the proof. 

We prove this claim by contradiction. Assume that no such c exists. 
Then there exists a sequence {u m , A m , a m , ZQ tTn , z m } such that for all m, f m 
attains its maximum in B am (zo tm ) at z m , and 

f{zm) := (cr m - dist(z m , z , m )) * u m (z m ) > m. 

Since (a m — dist(z m , zo, m )) 2 < ffm 2 < C for all m, we see that neces- 
sarily u m (z m ) -4 00. 

Let z be a system of Riemann normal coordinates centered at z m , so that 
the corresponding metric coefficients satisfy (g c )ij = &ij + 0(M 2 )- (As m 
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varies, these coordinate systems also vary, but there is no reason to include 

2 

this in the notation.) Set A m = (u m (z m )) n - 2 i ; we consider the dilated 
coordinate system £ = X m z, the corresponding sequence of metrics g m , 
where 

n 

9c = A~ 2 ^2 (.9c)ij((/^m)dCd( J := \~ 2 g m , 
and finally the dilated family of solutions 

n — 2-y 



By construction, v m (0) = 1 for all m, and 



A, 



t n — 27 n — 27 ~ 

5"! : — ^m 3c — Vm (Jm- 

We show below that g m and t> m are defined on an expanding sequence of 
balls on W n , and it is then clear that g m converges to the Euclidean metric 
uniformly in C°° on any compact subset. 

Let r m = \{a m - p(z m )), or equivalently, p m (z m ) + 2r m = a m . Then 

Cm P\%m) — Pm(.Zm) — Tm 

on the ball dist 9c (z, z m ) < r m , and hence on this same ball, 

n+27 

I \ / ( a m ~ Pm{Zm) \ "~ 2 '' , s . , s n ^2. 

U m {z) < -— Mm(Zm) < CM m (z m ), C=2 2 . 

V On. - Pm{z) J 

2 

The corresponding ball in rescaled coordinates contains : \Q\ < m n - 2 ~i}, 
hence has radius tending to infinity. By construction, v m {z) < c on this 
entire ball, and v m (0) = 1. Since these functions are uniformly bounded 
and satisfy the converging set of elliptic pseudo differential equations 

n + 2-, 

P° m v m = vZT 2 \ (4.15) 

we conclude using the local regularity theory (which is straightforward since 
v rn is bounded) that v m is bounded in C 2,a of every compact set, and hence 
we can extract a convergent subsequence. We thus obtain a smooth solution 
v to the 'flat' equation 

n+27 

— (A-gn^v = i;»- 2 7 inM n . (4.16) 
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Since each v m > 0, we see that v > 0, but v ^ since t> (0) = 1. There is a 
maximum principle for this equation \12\ Corollary 3.6] when < 7 < 1, so 
we conclude that v > on all of M. n . 

There is a complete characterization of positive solutions of (|4.16p . [T2"l 

§5]). They are the extremal functions for the embedding H^(M. n ) L n - 2 i (M. n ), 
and are necessarily of the form 



for some /j,, c > and zo € M n (these are well known "bubbles"). 

The argument is completed using Theorem 14.11 below, which states that 
a small ball in O m = S n \A m must have a concave boundary with respect to 
to g m for m sufficiently large. This is a contradiction to the already known 
limiting form of the v m . The proof of Proposition 11.31 is thus completed. 



Note that the previous arguments do not require that u is a weak solution 
in the whole S n . The only place where this strong hypothesis is required is 
in the following convexity claim: 

Theorem 4.1. In the same hypothesis as in Proposition \1.3l any open ball 
B (with respect to g c ) with B C £1, has boundary dB which is geodesically 
convex with respect to g. 

This result was proved in [21] for constant scalar curvature metrics, and 
also in the case 7 G (l,n/2) for locally conformally flat manifolds satisfy- 
ing some extra conditions by [22j . The crucial step is the application of the 
Alexandroff moving plane method. As we show here, the same ideas work in 
the fractional case. The moving plane method has been successfully applied 
to fractional order operators in [22] and [5], at least when the equation is 
rewritten as an integral equation. However, the proof in the present set- 
ing is simpler because of the equivalent formulation (|4.14p and the precise 
asymptotics (I4.19p . so we include the details for the reader's convenience. 
Our proof follows the classical arguments for the Laplacian by Gidas-Ni- 
Nirenberg in [TQ], [9]. 

\dx\ 2 

For simplicity, we denote P 7 := P\ . Let v be a distributional solution 




□ 



of 



P y v = v n - 2 i in M ; 



it 



(4.17) 
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We will apply the Alexandroff reflection with respect to the planes S\ := 
{x G W 1 : x n = A}. Let H\ := {x G W 1 : x n > A} be the hyperplane lying 
above S\. Given x = (x l , . . . ,x n ) € T,\, define x x to be the reflection of x 
with respect to the hyperplane S\, i.e., x x := (x 1 , . . . ,rr n ~ 1 ,2A — x n ). We 
also define v\{x) := v(x x ) and 

w\(x) := v x (x) - v(x). 

Note that the equation satisfied by v\ is the same as the satisfied by v. 
Although this fact is not clear for non-local operators, it is easily seen to be 
true in the Caffarelli-Silvestre extension (|4.14p . Then, by linearity, 



P lWx = v ^ - V n^ ) W eakly. (4.18) 

We will need a couple of preliminary results: 
Lemma 4.2. Let v be any function with asymptotics 

v(x) = |x| 27 " n (^ + ^2^ + 0(\x\ 2 )^j when \x\ -> oo, (4.19) 

for some a > 0. Then there exists Ao > such that for all A > Ao, 

w\(x) > for all x G T,\. 
Proof. This is just Lemma 2.2. in [10], and it does not use (|4,17p . □ 

Lemma 4.3. Let v a weak solution of (I4.17D , If for some A < Ao we have 
that w\{x) > but w\ ^ in Y,\, then 

w\(x) > in Sa and d n v{x) < on S\. (4.20) 

Proof. When v solves the constant scalar curvature equation, this is just 
Lemma 2.2. and Lemma 4.3 in [9]. In our case, we need a strong maximum 
principle and Hopf's lemma for the operator P 7 (see |12j . pQ). We know from 
(I4.18P that P^wx > 0. Since w\ > in T,\ (and is not identically zero), and 
since w\ vanishes on the boundary S\, the strong maximum principle gives 
that w\ > in all of T,\. On the other hand, Hopf's lemma implies that 
d n u)\ > on S\. Then, d n w\ = d n v\ — d n v = —2d n v, so we immediately 
have d n v < along S\. 

□ 
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Proof of Theorem \4-l\ Let g be a complete metric of constant positive 

4 

scalar Q 7 curvature on Q, C S n of the form g = u n ~ 2 i g c , and B an open 
ball in f2. Let S = dB be the boundary sphere. Fix any point p € S. 
Use stereographic projection to map O into f2 C M n so that p is mapped to 
infinity. Then S is transformed to a hyperplane 5, and the projected <9f2 
lies on one side of S, say below. Use linear coordinates (x , . . . E M n 
with 5 = {x n = 0}. 

By stereographic projection, the metric g transforms to a conformally 
fiat metric on M n , ^ = v n ~ 2 ~t \dx\ . Since the scalar curvature equation is 
conformally covariant, we also have 

n+2~f 

on C W 1 . Note that v is a weak solution of this equation on all of W l . 

Since the function u is smooth and strictly positive at p, the function 
v is regular at infinity, i.e. has the asymptotics ()4.19j) for some a > as 
\x\ — > oo. 

Step 1. Starting the reflection: Thanks to Lemma 14.21 we can initiate 
the reflection argument when A is sufficiently large. Note that the equation 
satisfied by v is not needed here since we have the precise behavior (|4.19p . 

Step 2. Continuation: We now move the plane S\, so long as it does not 
touch the singular set. Suppose that at some Ai > we have Wx x (x) > for 
all x € E^j , but w\ 1 ^ in . Then the plane can be moved further; more 
precisely, there exists some e > not depending on Ai such that w\ > in 
Y,\ for all A € [Ai — e, Ai]. 

We observe first that because of Lemma 14.31 we must have 

d n v < on S Al . (4.21) 

Next, the proof of our claim follows as in Lemma 2.3. in [10] by contra- 
diction. Thus, assume that there is a sequence Xj — > X\ and a sequence of 
points {xj}, xj E Y<\. such that 

w Xj (xj)<0. (4.22) 

Either a subsequence, which we again call {xj}, converges to Xoo G T*\ 1 or 
else Xj — > oo. In the first case, because of (|4.22p we must have d n v{x 00 ) < 0, 
thus contradicting (I4.2ip . So Xj — > oo. But in this second case may use the 
asymptotics for v from (|4.19j) . that imply 

\x ■ \ n 
Aj Xj 
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This is a contradiction to (|4.22p . 

Finally, note that in this process we never have w\ = since the existence 
of the singularity of v implies that it has no plane of symmetry. Hence the 
moving plane can be moved all the way to A = 0. 

Step 3. Conclusions: We have shown that T,\ can be moved to A = 0, 
and then w\(x) > for all x £ T,q and 

d n v(x) < for all x € So- 

4 

Since g = v n - 2 ~< \dx\ , the second fundamental form of any plane S\, A > 0, 
with respect to g is given by 

( ^-i^Lj] 

\ n — 27 dx n J 

The sign of d n v therefore implies that S\ is locally geodesically convex for all 
A > 0. When transferred back to the sphere, this shows that any round ball 
contained in has locally geodesically convex boundary, as claimed. □ 
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